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Abstract. Using the coupled-cluster method for infinite lattices and the exact diagonalization 
method for finite lattices, we study the influence of an exchange anisotropy A on the ground- 
£h state phase diagram of the spin-1/2 frustrated J\-Ji XXZ antiferromagnet on the square lattice. 

I 1 We find that increasing A > 1 (i.e. an Ising type easy-axis anisotropy) as well as decreasing 

A < 1 (i.e. an XY type easy-plane anisotropy) both lead to a monotonic shrinking of the 
parameter region of the magnetically disordered quantum phase. Finally, at A c « 1.9 this 
quantum phase disappears, whereas in the pure XY limit (A = 0) there is still a narrow region 
around J2 = 0.5Ji where the quantum paramagnetic ground-state phase exists. 
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A canonical model to study the interplay between frustration and quantum fluctuations in 
magnetic systems is the spin-1/2 Heisenberg antiferromagnet on the square lattice with nearest- 
. neighbour (NN) coupling Ji and frustrating next-nearest-neighbour (NNN) coupling Ji (J1-J2 

model), see, e.g., Refs. [1-8]. The recent syntheses of magnetic materials that can be well 
■ described by the J1-J2 model [9, 10] has stimulated further interest in the model. For the 

isotropic spin-1/2 J1-J2 model there are two magnetically ordered ground state (GS) phases 
at small and at large J2 separated by an intermediate quantum paramagnetic phase (QPP) 
without magnetic long-range order (LRO) in the region < J2 < J 2 2 , where J^ 1 ~ 0.4Ji 
and J 2 2 « 0.6 Ji. The GS at J 2 < J 2 Cl exhibits Neel LRO. The twofold degenerate GS at 
J2 > J 2 2 shows so-called collinear magnetic LRO. These two collinear GS's are characterized by 
a parallel spin orientation of nearest neighbours in the vertical direction and an antiparallel spin 
orientation of nearest neighbours in the horizontal direction [collinear-columnar (CC) state], and 
vice versa (collinear-row state). The nature of the transition between the Neel phase and the 
QPP as well as the properties of the QPP are still under debate [5-8]. 

Several extensions of the J1-J2 model have been studied recently, see, e.g., Refs. [11-21]. For 
instance, it was found that by increasing the spatial dimensionality from D = 2 to D = 3 the 
intermediate QPP disappears [12-14]. With respect to experimental realizations of the J1-J2 
model an exchange anisotropy could be relevant. Surprisingly, only a few papers so far have 
considered this issue. Some recent papers have discussed the special cases where (i) only the NN 
coupling Ji is anisotropic [20] or, alternatively, where (ii) only the NNN coupling J2 becomes 
anisotropic [19]. In materials it seems to be more likely that both couplings, Ji and J2, are 




Figure 1. The GS energy per spin scaled by its value for J2 = 0, E{J2)/E{J2 = 0), 
for anisotropics A = (XY), A = 1 (isotropic Heisenberg), A = 2 (Ising-type). (a) 
CCM: The LSUBn results with n = {4, 6, 8, 10} are extrapolated to n — ► 00 using E{n) = 
a + ai(l/n) 2 + a 2 (l/n) 4 . (b) ED: N = 36. 



anisotropic. The corresponding model is the square-lattice spin-^ J\-J% XXZ model 

H = Jx + 4s) + A S f S |) + J 2 £ + s?sl + Asfst) , (1) 

(iJ) 

where the first sum runs over all NN and the second sum runs over all NNN pairs. To our best 
knowledge the only study by other authors of such an anisotropic J1—J2 model with the same 
anisotropy in the J\ and the J2 term has be performed by Benyoussef et al. [11] using linear 
spin-wave theory. Moreover these authors considered A > 1 only. However, from early studies 
of Igarashi [2] of the isotropic J\-Ji model it is known that higher orders in the 1/s expansion 
become large near the quantum critical point (QCP), and hence results from the lowest order 
(i.e. linear) spin- wave theory become unreliable. As in our previous work [8,14,18,21] on J\-Ji 
models on the square lattice, we use here the coupled cluster method (CCM) complemented by 
exact diagonalisation (ED) of a finite square lattice of N = 36 = 6 x 6 sites (imposing periodic 
boundary conditions) to investigate the effect of exchange anisotropy. The CCM is an effective 
tool for studying highly frustrated quantum magnets [8,14,18,21-23], where, e.g., the quantum 
Monte Carlo method is not applicable due to the minus-sign problem. 

For the CCM treatment of the model (pQ), we use the classical GS (Neel at small J2 
and CC at large J2) as reference state |3>). Starting from these reference states the CCM 
employs the exponential parametrization = e 5 |<3?) of the quantum GS \^) where the 
correlation operator S contains all possible multi-spin-flip correlations present in the true GS. 
Naturally, S has to be approximated. We use the well-elaborated CCM-LSUBn approximation 
[8, 14, 18, 21-23] to calculate the GS energy per spin E and the sublattice magnetization 
per spin M. Since the LSUBn approximation becomes exact for n — ► 00, it is useful to 
extrapolate the 'raw' LSUBn data to n — > 00. There are well-tested extrapolation formulas, 
namely E(n) = ao + ai(l/n) 2 + a2(l/n) 4 for the GS energy per spin [14, 18, 21, 23, 24] and 
M(n) = b + biil/n) 1 / 2 + 6 2 (l/^) 3/2 for the sublattice magnetization [8,18,21]. We will not 
present more details of the CCM, but rather refer, e.g., to Refs. [8,14,22-24]. 

We start with the GS energy plotted in Fig. Q] for three characterictic values of the anisotropy 
parameter A. As mentioned above we use for the CCM calculations the Neel reference state at 
small J2, but the CC reference state at large J2. Hence, the CCM curves typically consists of 
two parts belonging to Neel and CC reference states. The curves for A = 1 shown in Fig. [TJa) 
illustrate clearly that there is a parameter region around J2 = 0.5 Ji where the CCM equations 
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Figure 2. Magnetic order parameter scaled by its value for J 2 = for anisotropies A = 
(XY), A = 1 (isotropic Heisenberg), A = 2 (Ising-type). (a) CCM: Sublattice magnetization 
M(J 2 )/M(J 2 = 0). The LSUBra results for M with n = {4,6,8,10} are extrapolated 
to n — > 00 using M(n) = bo + &i(l/n) 1//2 + 62(l/n) 3 / 2 . The inset shows the parameter 
region around J 2 = 0.5Jj with an enlarged scale, (b) ED: Square of the order parameter 
(m+) 2 (J 2 )/(m + ) 2 (J 2 = 0) for N = 36. 



using classical Neel and CC reference states do not have real solutions. As a consequence, the 
corresponding pair of GS energy curves for the Neel and CC phases do not intersect one another. 
This behaviour yields preliminary evidence for the opening up of an intermediate phase between 
the Neel and CC phases. By contrast, for A = and A = 2 the corresponding pairs of GS 
energy curves almost cross one another giving preliminary evidence that for strong anisotropy 
a direct first-order transition between both semiclassical magnetic phases may occur. From the 
ED data it is also evident, that the behaviour of the GS energy for the isotropic model, i.e. 
A = 1, differs from that for A = and A = 2. 

Next we present in Fig. [2] the sublattice magnetization M calculated by the CCM for N — ► 00 
and the square of an order parameter defined as (m + ) 2 = ^Zij I ( s i ' s j) I [25] calculated by ED 
for iV = 36. More data for M can be found in Ref. [21]. While M is finite in the magnetically 
ordered phases but vanishes in the intermediate QPP, we have always finite values for (m + ) 2 
for finite N = 36. We use the CCM results for M to define the QCP's (J 2 C \ A Cl ) and (J 2 C2 , A C2 ) 
as that points where the magnetic order parameter [Neel at (J^A^) and CC at (J 2 2 , A C2 )] 
becomes zero. From Fig. [2] it is obvious that the intermediate QPP is largest for A = 1 (the 
CCM estimates for J 2 Cl and J 2 C2 are J 2 Cl « 0.44 . . . 0.45 Ji and J 2 2 « 0.58 . . . 0.59 Ji for A = 1, cf. 
Refs. [8,21]). Both types of anisotropy lead to a stabilization of magnetic LRO. The ED data for 
(m + ) 2 support these findings. From Fig. [^b) it can be seen that the parameter region of small 
values of (m + ) 2 around J 2 = 0.5Ji is significantly broader for A = 1 than for A = 2 and A = 0. 
From the inset of Fig. [^a) it is also obvious that the pair of CCM order-parameter curves (Neel 
and CC) for A = 2 intersect one another at a value M > 0. This behaviour can be interpreted 
as an indication of a first-order transition between the magnetically ordered Neel and CC phases 
located at the crossing point, for a detailed discussion of this issue see also Refs. [14,21]. 

The influence of the anisotropy on the correlator (so • sr) is illustrated in Fig. Eta) for 
J 2 = 0.45Ji, i.e. near the QCP J 2 X where the Neel LRO breaks down for A = 1. We see that 
(so • sr) decays most rapidly for A = 1. For the largest separation R m &x = V^8 present for 
iV = 36 sites, the correlator (so • s^ max ) for J 2 = 0.45Ji is reduced by frustration by a factor of 
0.25 for A = 1, whereas the corresponding factor is only 0.52 (0.40) for A = 2 (A = 0). 

Finally, we summarize our findings on the GS magnetic ordering in Fig. EKb) where the 
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Figure 3. (a): ED results for the spin correlation function scaled by its value for J<i = 0, 
(so • s.r)(J2)/(so • s#)(J2 = 0) versus separation R for A = 0, 1, and 2 and J2 = 0.45 Ji. (b) GS 
phase diagram of the spin- 1/2 J\-J% XXZ model on the square lattice calculated by the CCM. 



GS phase diagram is shown. The solid lines in Fig. E^b) represent those points in the (A, J2) 
parameter space where the Neel and the CC order parameter calculated by the CCM vanish. 
Increasing the anisotropy leads to a monotonic shrinking of the region of the QPP. For the easy- 
axis anisotropy all three phases meet at a quantum triple point at (A c ~ 1.9, J| ~ 0.52), i.e. 
the QPP disappears completely for A > 1.9. Similarly, for the case of the easy-plane anisotropy 
a second quantum triple point occurs at (A c ~ — 0.1, J| ~ 0.50). Outside the area hemmed 
by the solid lines there is a direct first-order transition between the Neel and the CC phase, as 
indicated by the dotted lines. 
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